Quantum signatures of self-trapping transition in attractive lattice bosons 
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We consider the Bose-Hubbard model describing attractive bosonic particles hopping across the 
sites of a translation-invariant lattice, and compare the relevant ground-state properties with those 
of the corresponding symmetry-breaking semiclassical nonlinear theory. The introduction of a suit- 
able measure allows us to highlight many correspondences between the nonlinear theory and the 
inherently linear quantum theory, characterized by the well-known self-trapping phenomenon. In 
particular we demonstrate that the localization properties and bifurcation pattern of the semi- 
classical ground-state can be clearly recognized at the quantum level. Our analysis highlights a 
finite-number effect. 



I. INTRODUCTION 

Over the last few years the Bose-Hubbard (BH) model 
[l| has attracted considerable attention, owing to its re- 
alization in terms of ultracold atoms trapped in optical 
lattices 0, Among other aspects, the interest in the 
relation between the BH model with attractive interac- 
tions and the nonlinear theory that can be considered its 
semiclassical counterpart [4-6] was rekindled One 
of the most striking issues in this respect lies in the sym- 
metry properties of the typical states of the theory. The 
translation symmetry of the quantum states appears to 
be at odds with the spontaneous symmetry breaking of 
the localized, self-trapped semiclassical states, although 
the same regimes can be identified in the two theories 
0, Hj]. These apparently conflicting features have been 
reconciled in a recent work [Toj| where it is discussed how 
a single measurement in a (inherently linear) quantum 
system can give rise to the localized, symmetry-breaking 
result typical of the corresponding nonlinear theory. The 
point in question is illustrated in some detail for the well 
known case of a lattice comprising only two sites, and 
a few further examples are given for a larger lattice, at 
a single value of the effective parameter governing the 
system. 

Here we extend such analysis in a twofold way. On the 
one hand, we introduce a suitable notion for the width of 
a localized state, thereby providing a more quantitative 
tool for comparing the quantum measurements with the 
relevant semiclassical predictions. On the other hand, 
unlike Ref. [To|, we do not limit our simulations to the 
localized regime, but we perform a systematic compari- 
son, exploring a range of effective parameters including 
the semiclassical bifurcation of the system ground state 
pH fl2j . The analysis of the localization width shows that 
the semiclassical derealization transition corresponds to 
a crossover in the quantum system, which highlights a 
finite-population effect. The sharp transition is recovered 
only in the limit of infinite boson filling, where the clas- 
sical results are exact. 



Also, the range of effective parameters and lattice sizes 
we explore includes those highlighting the change in the 
nature of the bifurcation occurring at the semiclassical 
level when passing from five-site to six-site lattices [l3|, 
Il4j ]. As we discuss in the following, a signature of this 
qualitative change is apparent at the quantum level. 

The layout of this paper is the following. In Sec- 
tion |TT] we review the second-quantized model for attrac- 
tive bosons hopping across the site of a lattice and the 
nonlinear theory representing its semiclassical counter- 
part. We discuss the localization/delocalization transi- 
tion occurring in the ground- state of the latter, revisit- 
ing the results discussed in Ref. [To[. Also, we provide 
a quantitative measure for localization, which is readily 
extended to the quantum theory. 

Section HIT1 contains our results, which are based on the 
numerical calculation of the ground-state of the quantum 
system. This is obtained by means of Lanczos diago- 
nalization and population quantum Monte Carlo simula- 
tions. We start with a detailed analysis of the specific 
case of a three-site lattice (trimer), for which the local- 
ization/delocalization transition can be appreciated by a 
direct visualization of the structure of the ground-state. 
We then make a systematic analysis of the localization 
width for several lattice sizes and boson fillings, high- 
lighting the agreement with the semiclassical result for 
large fillings. This also shows that the qualitative change 
in the bifurcation pattern of the classical theory taking 
place for lattices containing more than five sites is recov- 
ered at the quantum level. 

In Section [TV] we discuss the connection between the 
quantum and classical theory based on the representation 
of the quantum ground-state as a su(L) coherent state Q 
(also known as Hartree wave function Q). This allows 
us to employ the semiclassical ground state to construct 
an approximation to the quantum ground which includes 
some of the finite population effects mentioned above. 
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II. THE SYSTEM 



We consider a simple Bose-Hubbard model 
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with attractive interactions, U, J > 0. The operators d/- 
and (ij create and destroy bosons at lattice site j, respec- 
tively. The lattice comprises L sites and is periodic, so 
that site label j = L is to be identified with j = 0. The 
(positive) parameters U and J account for the relative 
strength of the interaction and kinetic term, respectively. 

The BH Hamiltonian clearly commutes with the to- 
tal number of bosons in the system TV = . hj , with 

hj — a}-aj. This means that the eigenstates of Eq. (pp) 
are characterized by a well defined total number of bosons 
N, i.e. that H can be studied within a fixed-number 
subspace of the infinite Hilbert space. The size of this 
fixed-number subspace is finite, s = ( iV+ ^ _1 ), but it 
becomes computationally prohibitive already for modest 
lattice sizes and boson populations. 

Hamiltonian (pp) is clearly translation invariant. By 
virtue of the Bloch theorem, its eigenstates — and in 
particular its ground state — are delocalized over the 
entire lattice. This means nj = (^f\hj\^f) = (^\h £ \^} = 
rig for any j and t. 

The semiclassical nonlinear theory corresponding to 
the BH model can be obtained by changing the site op- 
erators into C-numbers, a}- i/jj, whose square modulus 
has a natural interpretation as the local boson popula- 
tion, hj ^> IV^'I 2 - The semiclassical Hamiltonian thus 
modeled onto Eq. (pp) results in a dynamics governed by 
the so-called discrete nonlinear Schrodinger equations or 
discrete self-trapping (DST) equations, 

ij>j = -U\^j\ 2 ^j - J(^+i (2) 
where we set H = 1 . The stationary solutions ipj (t) = 



obey the fixed point equations 
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Similar to the quantum case, the semiclassical Hamil- 
tonian, as well as equations of motion (|2]) and (|3]), are 
translation invariant. However, their nonlinear nature al- 
lows of symmetry-breaking solutions. In particular, it is 
easy to check that Eq. ([3]) always has a uniform solution, 
4>j — y/N/L, independent of the value of the parameters 
U and J. Its frequency and energy are lj = —UN/L — J 
and E = —UN/2L, respectively, the latter being the low- 
est possible energy only for sufficiently large values of the 
effective parameter 



J 
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If, conversely, the effective interaction energy UN pre- 
vails over the hopping amplitude J, so that r < ra, 



the lowest-energy solution of Eq. (j3j) is localized about 
one lattice site, thus breaking the translation invariance. 
This phenomenon is referred to as (discrete) self-trapping 
In general, a second critical value 
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can be recognized for the uniform solution of Eq. (j3j) such 
that for r > r s the uniform solution is dynamically stable, 
whereas it is unstable for r < r s [12|. 

For sufficiently large lattice sizes, L > 6, the equality 
applies in Eq. ([5]). That is, the uniform solution becomes 
simultaneously stable and the ground state of the system 
as soon as r > r s = t&. On smaller lattices the strict 
inequality applies, r s < ra, and the low-energy solutions 
of exhibit a more complex bifurcation pattern [l3l-[l5| . 

A signature of the localization occurring at the semi- 
classical level for large values of the interaction can be 
recognized also at the quantum level. The quantum 
ground-state can be seen a symmetric superposition of 
L states, each localized at one of the L sites of the lattice 
and closely resembling the symmetry-breaking semiclas- 
sical ground-state. At the boson fillings N/L correspond- 
ing to Hamiltonians that can be analyzed by means of 
Lanczos diagonalization, this resemblance is strong only 
for very small values of the effective parameter r, while 
it is washed out by quantum fluctuations at larger r's [§[. 

In Ref . [HJ the connection between the uniform quan- 
tum ground state and its localized semiclassical coun- 
terpart is further discussed. There it is remarked that 
the experimental measurement of the observable n = 
(ni, n2, • • • , ul) is likely to produce an outcome in strong 
agreement with the semiclassical result. Indeed such a 
measurement selects a Fock state of the direct space, i.e. 
an eigenstate of n, 



fvj\ 



|0», 
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where |0)) is the vacuum of the theory, i.e. a>j\0}} for 
all j = 1, 2, . . . , L. Each Fock state is selected with a 
probability distribution given by P(y) = |c^| 2 , where the 
c^'s are the coefficients in the expansion 



I*) = £W) 



(7) 



The prime on the summation symbol signals that the 
Fock states involved in the expansion belong to the same 
fixed- number subspace, J2<i v j = ^« 

The Authors of Ref. [10] first of all recall that for 
a two-site lattice, L = 2, the probabilit y d istribution 
P(y) is double peaked when r < Td [HI, [I?} , and ob- 
serve that the mirror-symmetric Fock states correspond- 
ing to the peak probabilities reproduce the populations 
of the symmetry-breaking semiclassical ground states: 

KkO = (IV'ilM^I 2 ), iyi,»2) = (hfc| 2 ,hM 2 ). Also, 
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by exactly diagonalizing Hamiltonian (pQ) they show that 
the width of the probability peaks decrease with increas- 
ing total number of bosons. 

The Authors' observation is further illustrated by con- 
sidering the case of N = 256 bosons on a lattice compris- 
ing L = 16 sites. A few Fock states \v) are sampled from 
a quantum Monte Carlo calculation and, after suitable 
translations, graphically compared with the semiclassical 
result for the local populations | 2 , IV^I 2 ? • • • , l^z, | 2 )- 
The importance-sampled quantum configurations agree 
reasonably well with the semiclassical result, although 
quantum fluctuations are clearly recognizable. 

In summary, Ref. [loj demonstrates how a measure- 
ment in an inherently linear quantum theory could pro- 
duce a typical outcome of a nonlinear theory, i.e. a lo- 
calized soliton-like ground-state. 

In the following we further investigate the connection 
between the ground-state of Eqs. (pQ) and (|3]), making the 
comparison more quantitative with the aid of a suitable 
observable. Also we explore a range of effective parame- 
ters r including the localization/delocalization threshold 
Td, and discuss the signature of this semiclassical tran- 
sition at the quantum level. In particular we illustrate 
that the difference in the semiclassical bifurcation pat- 
tern distinguishing small lattices from those comprising 
more than six sites [l3|, [H| is apparent also in the quan- 
tum data. 



A. Soliton width 

As we recall above, for r < the ground-state of the 
DST equations ([3]) becomes localized at one of the lat- 
tice sites, assuming a soliton-like density profile. Since 
the soliton peak can be localized at any of the L lattice 
sites, the ground-state is L-fold degenerate. This spon- 
taneous symmetry breaking entails from the nonlinear 
nature of Eq. (j3]). After recalling that we are consider- 
ing a cyclic lattice, i.e. periodic boundary conditions, we 
can estimate the (square) width of a localized solution of 
Eq. (j3j) as 
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where Xj = cos ^-j and z/j 



sin are the coordinates 



of the j-th site of the ring lattice [18|, rij 
boson population at that site and 



|^-| 2 is the 
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are the coordinates of the center of mass of the boson 
distribution. When the semiclassical solution is uniform, 
rij = \i/jj\ 2 = N/L, the center of mass is at the center 
of the ring lattice, x cm = y cm = and the width attains 
the maximum possible value, w = 1. At the opposite 




FIG. 1: Width of the semiclassical ground-state according to 
Eq. ((8]) for some lattice sizes, from left to right 3 < L < 
7. Notice that the transition from a localized (w < 1) to a 
uniform state (w = 1) is discontinuous for L < 6 [13, [3] • The 
dotted lines at the critical values ra are guides to the eye. 



limit, the center of mass coincides with the lattice site 
at which the entire boson population is confined. In this 
situation the width attains its minimum value, w = 0. 
The width of the lowest solution to Eq. ((3]) is plotted in 
Fig. [I] for some small lattices, 3 < L < 7. Note that the 
different bifurcation pattern characterizing the smaller 
lattices [l3|, [H|, L < 6, is mirrored in the discontinues 
character of the width of the ground state. In general 
an increase of r < Td results in an increase of the width 
of the localized state. However, the derealization is at- 
tained continuously at r = Td only for L > 6. For smaller 
lattice sizes the width of the localized ground state has 
an upper bound smaller than 1, and derealization is at- 
tained catastrophically as r > Td- 

As we discuss above, a measurement of the quan- 
tum observable n selects a Fock state, i.e. a set of in- 
teger occupation numbers Vj = (fij), with probability 
P(v) = \ca\ defined by Eq. Eq. (j8j) lends itself to 
the estimate of the width of the selected Fock state as 
well, provided that Vj is plugged in n^, instead of |^j| 2 - 



Taking into account that x\ 



y 2 



1 and that the occu- 



pation numbers in the Fock states add up to A 7 ", after a 
few manipulations Eq. (|8]) becomes 
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Therefore, the average of a large number of measure- 
ments of w tends to the quantum observable 



= 1 



v 

1 ^ 



N 2 

1 

2 



cos 



L 



(hjhi) 



(ii) 



where 

^) = ^E e " (H) ^- ( 12 ) 

We note that the quantity in Eq. (fl2j) is formally similar 
to the static structure factor for a ID linear lattice pjl 

III. RESULTS 

In the following we show that the average width of a 
quantum ground state, as estimated by Eq. ([TT]h repro- 
duces the semiclassical results shown in Fig. [TJ the agree- 
ment improving as the boson population in the lattice is 
increased. We thereby make the analysis of Ref. [Toj| 
more systematic and quantitative. Furthermore, we do 
not limit our investigation to the region where the semi- 
classical solution is localized, but explore the transition 
region as well. We do this for relatively small lattices, 
but consider one example, the six-site lattice, for which 
the bifurcation pattern at the transition is qualitatively 
similar to larger lattices [l3|, Q • 

In order to evaluate the quantity in Eq. ([TT]) we need to 
obtain the ground-state of Hamiltonian (pQ). We achieve 
this by means of two different numerical approaches. 
When the size of the relevant Hilbert space is sufficiently 
small, we employ a Lanczos diagonalization algorithm. 
Otherwise, we resort to a stochastic method, the so-called 
population quantum Monte Carlo algorithm (2lj |. In both 
cases we exploit the symmetry granted by the commuta- 
tor [H, N] by working in the canonical ensemble, i.e. by 
considering only the occupation- number Fock states \v) 
relevant to a given total number of bosons, TV = J2j v r 
In order to reduce further the size of the matrix to be an- 
alyzed by the Lanczos algorithm, we also take advantage 
of the translation symmetry of the system. We first of 
all gather the fixed-number Fock states into equivalence 
classes determined by lattice translations, which allow to 
define the reduced basis: 

|,;) = — ^M^l.l (is) 
VMa j=1 

In the r.h.s. of this equation \v) represents any member 
of an equivalence class, which determines all of the other 
Ma — 1 members of the same class through the displace- 
ment operator!) such that DajD~ l = The number 
Ma of states in a class is in general a divisor of the lattice 
size L. The size of the reduced basis comprising states 
of the form ([T3]) is the same as the size of the subspace 
formed by the quasimomentum Fock states \q) such that 
J2k=i^^k = k L. with k G Z i.e. roughly a factor L 
smaller than the fixed-number Fock space. 

Before discussing the results for the average quantum 
width, Eq. (fTT]h we analyze the ground state of three-site 
lattice in some detail. 




FIG. 2: Representation of the probability distribution P(y). 
Left: the dark lines demonstrate the (flat) surface of the 
occupation- number Fock space relevant to a trimer containing 
N bosons. Center: the same surface as viewn from direction 
(001), i.e. projected onto the (^1,^2) plane. Right: the same 
surface as in the rightmost panel as viewn from the direction 
normal to it, (111). Notice that this view highlights the three- 
fold symmetry of the domain of the distribution probability. 



A. The trimer 

For lattices comprising just three sites the probability 
distribution P{y) for the Fock state \v) = ^1,^2,^3) se- 
lected in an experiment measuring n = (ni, 77,2,^3) can 
be conveniently represented as a two dimensional density 
plot. This is made possible by the total number conser- 
vation, which constrains one of the occupation numbers, 
say vs , to a value depending linearly on the two remaining 
occupation values ^3 = N — v\ — v^. One can therefore 
regard the probability distribution as a function of the 
latter alone, P{y\, v<i, N — v\ — v 2 ). 

The portion of the occupation-number Fock space rel- 
evant to a trimer containing N boson is illustrated in the 
leftmost panel of Figure [2j The same surface as seen from 
two different points of view is shown in the remaining two 
panels. In the central panel the point of view is on the 
direction (0,0, 1), which results in a projection onto the 
(^1,^2) plane. In the leftmost panel the point of view 
is on the direction (1,1,1) normal to the surface under 
investigation. This highlights the three-fold symmetry of 
the surface, which is not apparent in the previous view. 
Since we are interested in the symmetry of the system, 
we adopt the second point of view when representing the 
probability distribution P(v). 

Figure [3] shows the density plots of P(y) for several 
values of the effective parameter r defined in Eq. (J2J). 
The representation is the same as in the rightmost panel 
of Fig. [2j but in order to compare different total popula- 
tions, the tick labels refer to Vj/N instead of Vj. The left 
and right columns correspond to two different choices for 
the total population, N = 300 and N = 900, respectively. 
Note that the probability density undergoes a qualitative 
change as r crosses the region of the semiclassical der- 
ealization point, Td = 0.25. Below this threshold P{v) 
features three symmetrically positioned peaks. Notice 
that the probability peaks occur at the same positions 
as the three equivalent semiclassical ground states. The 
white cross symbol signals the position of one of such 
semiclassical ground- states, namely the one with v\ = v 2 
and us = N — 2 v\ . Owing to its symmetry, the three- 
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modal distribution always produces a symmetric expecta- 
tion value of the occupation numbers, (hj) = N/3. How- 
ever, a single measurement of (ni, ^2,^3) is extremely 
likely to produce a symmetrically broken outcome very 
similar to the semiclassical result. 

As r increases above ra, the three- modal probability 
distribution becomes monomodal. In this situation the 
symmetric expectation value (hj) = N/3 is very similar 
to the outcome of a single measurement, in agreement 
with what happens at the semiclassical level. 

Notice that there is a very narrow range of effective 
parameters in which the probability is essentially four- 
modal, featuring a central peak surrounded by three sym- 
metrically positioned peaks. The value of r such that 
the central peak has the same height as the peripheral 
peaks could be used to define the quantum counterpart of 
the derealization threshold t&. This situation is demon- 
strated in the density plots in the fourth row of Fig. [H 
These plots refer to slightly different values of r, depend- 
ing on the different total populations. As the population 
increases, this quantum threshold approaches the semi- 
classical one, Td = 0.25. The remaining rows of Fig. [3] 
refer to the same values of r. 

The quantum threshold condition discussed above, 
max^P(z7) = y, y), requires a tomography of 

the quantum state. A substantially equivalent approach 
makes use of the average width defined in Eq. (fTTj) . 



B. Average width 

In this section we demonstrate how the measured aver- 
age width for a quantum system, Eq. ([TT]h reproduces the 
semiclassical results shown in Fig. HJ Once again we first 
of all consider the three-site lattice. Figure |4] shows the 
average quantum width, Eq. (fTT|) for a lattice compris- 
ing L = 3 sites and for different total boson populations, 
N = 300, 600 and 900 particles. The semiclassical result 
appearing in Fig. [T] is also plotted for purposes of com- 
parison. Notice first of all that the range of r's in Fig.|4]is 
much smaller than that in Fig. [TJ and outside this range 
the quantum results corresponding to the three consid- 
ered populations agree extremely well with the semiclas- 
sical findings. The differences between the three sets of 
quantum data and the semiclassical ones become per- 
ceptible in a narrow region surrounding the semiclassi- 
cal localization/delocalization threshold, signalled by the 
vertical dashed gray line. Unlike its semiclassical coun- 
terpart, the average quantum width is a continuous func- 
tion of the effective parameter r. However, its derivative 
becomes extremely large in a very narrow region in the 
proximity of the semiclassical threshold. The width of 
such region decreases with increasing boson population 
so that, in the limit N — >• 00 the semiclassical disconti- 
nuity corresponding to the catastrophe in the bifurcation 
pattern of the semiclassical ground state [§, [l3j becomes 
apparent. The points on the almost vertical stretch of 
the quantum curves correspond to four-modal probabil- 
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FIG. 3: Density plots of the probability distribution P{v) 
for some values of the effective parameter r, Eq. Q. Left: 
TV = 300. Right: N = 900. The probability distribution is 
represented as described in the rightmost panel of Fig. El but, 
in order to compare different total populations, the tick labels 
refer to fj/N instead of h>j. A white cross symbol signals the 
position of one of the three equivalent semiclassical ground- 
states (see text). The density plot is scaled with respect to 
Pa r = max,-? P( is) . 
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FIG. 4: Average quantum width Eq. (fTT|) for a trimer con- 
taining 300, 600 and 900 bosons. The semiclassical result is 
also plotted for purposes of comparison. 
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FIG. 5: Average quantum width Eq. ((11)) for a L — 5 lattice 
and several boson populations. The solid line black line for 
N = 60 was obtained by means of Lanczos diagonalization, 
while all the data points are population QMC results. The 
dotted lines are mere guides to the eye. The semiclassical 
result is also plotted for purposes of comparison. 



ity distributions, like in the fourth row of Fig. [2j 

The same qualitative behaviour is observed on lattices 
comprising four- and five-sites. The data corresponding 
to several boson populations are shown in Fig. [5] for the 
latter case. Even exploiting all of the system symme- 
tries, the Hamiltonian matrix for the L = 5 lattice be- 
comes extremely challenging for the Lanczos algorithm 
around fillings of the order of 12 particles per site. We 
have nevertheless been able to tackle larger fillings by re- 
sorting to a stochastic approach, the so-called population 
QMC algorithm (2l|. In Fig. [5] we also plot the results 
for N = 60 as provided by the Lanczos algorithm. The 
comparison with the data points for the same population 
demonstrates the reliability of our population QMC. 



FIG. 6: Average quantum width, Eq. (fTTj) . for a 6-site latice 
at different boson populations. The solid line corresponding 
to 30 (36?) bosons has been obtained by means of Lanczos 
exact diagonalization. The data points on top of it have been 
obtained by means of population QMC 21], as well as the data 
points for the larger populations. The dotted lines are mere 
guides to the eye. The semiclassical result is also plotted for 
purposes of comparison. 



As already observed in the trimer case, Fig. |4j outside 
a narrow range of r's the quantum data corresponding 
to sufficiently large boson fillings overlaps very well with 
the semiclassical result. An incipient discontinuity is rec- 
ognized in the proximity of the semiclassical threshold, 
which becomes more and more evident with increasing 
total population. The very same behaviour is oberved in 
the L = 4 case, not shown. 

The qualitative change at the transition taking place 
for lattice sizes larger than 5, evident in Fig.Q] is clearly 
recognizable also at the quantum level. In order to illus- 
trate this we consider the threshold situation, i.e. a lat- 
tice comprising L = 6 sites. The data points in Fig. [6] — 
also obtained by means of the population QMC algorithm 
- illustrate the behaviour of the quantum width at sev- 
eral boson populations. The solid green line once again 
demonstrates the agreement between QMC and Lanczos 
algorithm at filling N/L = 6. As the boson population 
is increased the data points for the quantum width get 
closer and closer to the semiclassical result which, at vari- 
ance with smaller lattices, is continuous at the transition 
point. 

Before concluding this paper we discuss in deeper 
detail the connection between the quantum theory of 
Hamiltonian (pQ) and its semiclassical counterpart. 



IV. su(L) COHERENT STATES 

The DST fixed-point equations ([2j) can be derived by 
assuming that the system is well described by a trial state 
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of the form 
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|0»- 



(14) 



This was first sketched in Ref. p], where the trial state 
(p~4|) was referred to as Hartree wave function, and subse- 
quently recast in terms of the time-dependendent vari- 
ational principle in Ref. [8], where Eq. (p3J) was rec- 
ognized as a su(L) coherent-state (CS). After extrem- 
izing a suitable functional, both approaches result into 
Eq. ([2j), except for a correction factor (N — l)/N cropping 
out in the interaction term. This comes about because 
Eq. (p3J) is an eigenstate of the total number operator, 
N\ij;) = N\ip), and is consistent with the expected ab- 
sence of interaction when only one boson is present in 
the system. This correction provides a better agreement 
between quantum and semiclassical theories for small bo- 
son populations but can be safely neglected at most 
of the fillings considered in the present paper. 

The probability that a measurement of n on the 
CS (fT4j) selects a Fock state |z7), Eq. (j6j), is given by 



P{V) = \{v\i,)f 



N\ 



n 



(15) 



(see e.g. Ref. [22j). Note that the above result applies 
when the total population of the Fock state is the same 
as the CS, (v\N\v) = Vj = TV. By using standard 
methods (see e.g. Ref. [23j) it is easy to appreciate that 
the above probability distribution is sharply peaked at 
the Fock state best reproducing the set of semiclassical 



occupation numbers, Vj w (^fi^ijj) = \^j\ 2 (see also the 
discussion about Fig. [7] below). Since the dynamical vari- 
ables ipj are determined by the DST equations (|2j), the 
CS breaks the translation symmetry below the localiza- 
tion threshold. That is, denoting D the matrix producing 
a ciclic shift of the vector entries, (Dv)j = Vj+i, one gets 
P(D k v) ^ P(v) for k^L. 

In Ref. [g| it was shown that, for sufficiently strong in- 
teractions, the above symmetry-breaking CS of is well de- 
scribed by an uniform superposition of the lowest-energy 
states of the L quasimomentum blocks of Hamiltonian 
(P). Here we take a somewhat complementary stand- 
point, and consider a uniform superposition of L equiva- 
lent, symmetry-breaking CS 



L L 



D 



(16) 



where the displacement quantum operator D is defined 
at the beginning of Sec. Mil The L (nonuniform) entries 
of the vector 4> are normalized such that {(4>\4>Jj = {4>\4>) = 



N -1 \(f>j\ 2 = 1, are yet to be determined. The optimal 
set of entries is such that the above symmetrized state 
attains the minimum energy. This entails 
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(17) 
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where A is a Lagrange multiplier enforcing the constraint 
on the norm. More explicitly this means 



ae^-w = -E 



N - 1 



TV- 1 
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N-2 tt i* i2 



U< t>hK+r + JU r 1 (0/i+r-l + <fe+r+l) 
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(19) 



and 

E n ^ = 1 ( 2 °) 

r 

where we introduced the shorthand notation 
II r = TV -1 J2 S (frs^s+r- Since nonuniform 0's 
yield IIo > IU, the only possible solution of Eq. ([20]) in 
the N — >• oo limit is IIo = 1, i-e- the usual normalization 
condition J2j l^il 2 = N. This entails that 11^ 7V ^° § r ^^ 
which considerably simplifies Eq. (fT9]l 

(A - S)(/> h = -^-±U\0 h \ 2 h - J(</> h+1 + (21) 



where £ = ($\H\$). Comparing Eqs. (|3]) and Eq. (I2T1) 
shows that, in the large population limit, the parameters 
cj) to be plugged into the symmetrized CS (fl6j) do not 
need to be obtained from the (numerically demanding) 
minimization inherent in Eqs. (fTTj) and (fl8]h but can be 
replaced by the (much more easily determined) normal 
modes of the dynamical DST equations ((2]), i.e. the so- 
lutions of Eq. (j3j). The comparison of Eqs. (j3j), (j4]) and 
(|2T|) shows that semiclassical parameter corresponding to 
Eq. (f2T|) must be rescaled by a factor jfzj- 

In Fig. [7] we show some density plots for the proba- 
bility distribution P(y) = |(^|</>))| 2 corresponding to the 




FIG. 7: Density plots of the probability distribution P(y) = 
|(z7|0))| 2 corresponding to the symmetrized coherent state in 
Eq. (|16[) for a trimer containing N = 300 bosons. We choose 
the same values of the effective parameter r, Eq. J4} as in the 
first, third and fourth panel from top in the left column of 
Fig. [3] A white cross symbol signals the position of one of the 
three equivalent semiclassical ground-states (see text). The 
density plot is scaled with respect to Pm = max^ P(y)> 



symmetrized coherent state in Eq. ([T6]) for a three-site 
lattice containing N = 300 bosons. The parameters </> 
are not obtained by solving Eqs. (fT7|) and (fT8|) . Instead, 
in keeping with the above discussion, we use the solu- 
tion of the DST fixed-point equation (J3j), given the fairly 
large boson filling. The parameter r controlling the rel- 
ative strength of the tunneling and interaction energy, 
Eq. (J2J), has been chosen as in the first, third and fourth 
panels in the left column of Fig. [3j respectively. As we 
mentioned earlier, the probability density ([T5]) for the 
su(AZ') coherent state ([T4]) is strongly peaked at the occu- 
pation numbers {vj} best reproducing the semiclassical 
local populations {|</>j| 2 }. It is then clear that, in the 
localized regime, the density of the symmetrized state 
(fT6|) is qualitatively similar to that of the quantum state 
analyzed in Fig. 03 It features three peaks, each rele- 
vant to one of the Fock states corresponding to the tree 
equivalent, symmetry-breaking solutions of Eq. (|3j). For 
sufficiently large interactions (small r's) the similarity is 
striking, as it was already pointed from a different per- 
spective in Ref. [§| . The differences between the quantum 
ground state in Fig. [3] and the symmetrized CS in Fig. [71 
become evident in the vicinity of the semiclassical thresh- 
old. In particular, at variance with the former, the latter 
is structurally unable to give rise to a four-modal distri- 
bution, as it is clear after comparing the central panel of 
Fig. and the third panel from top in the left column of 
Fig. El 

Above the derealization threshold there is no need of 
symmetrizing the CS in Eq. (fl4|h because the ground- 
state solution of Eq. (j3j) becomes translation invariant, 

= yj N/M, which makes the CS and the corresponding 
probability density r-independent. As it is clear from 
the comparison between the rightmost panel in Figs. 
and the bottom panel in the left column of Fig. [3j the 
similarity between the CS and the quantum ground state 
is already rather satisfactory at r = 1. It improves with 
increasing r, and becomes perfect in the noninter acting 
limit. 

We conclude by observing that the (quantum) width 
of CS ([T6]) as calculated with Eq. ([TT]) coincides with 
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FIG. 8: Comparison of the (squared) width of the semiclas- 
sical ground state, the coherent state and the exact quantum 
ground-state for a six-site lattice. The latter results refer to 
total population of N = 36 bosons, a Note that the three 
curves overlap in the small-r region, and that the coherent 
state provides a better approximation in the large-r limit at 
the relatively small filling considered. 

the semiclassical value Eq. (j8j) only in the large- TV limit. 
Straightforward calculations show indeed that 

«KI#) = ^™({l^| 2 }) ( 22 ) 

The population-dependent prefactor in the previous 
equation correctly allows of the expected derealization 
when one single boson is present in the system. Also, 
it provides a better agreement with the exact result in 
the large-r limit. As we discuss above, the same pref- 
actor affects the semiclassical parameter r to be used in 
Eq. (|2T]h defined in Eq. This correction makes it so 
that the width of the symmetrized CS is closer to that of 
the quantum ground state in the small-r (large interac- 
tion) regime. This should be clear from Fig. [5J and shows 
that the su(L) CS approach to the semiclassical theory is 
more effective in capturing the effects arising from finite 
size than the simple replacement of lattice operators by 
C-numbers used in deriving Eq. (|2]). 



V. CONCLUSIONS 

In this paper we perform a systematic analysis of 
ground-state properties of a system of attractive lat- 
tice bosons, highlighting the correspondences between 
the (inherently linear) quantum theory for this system 
and its nonlinear semiclassical counterpart. Our anal- 
ysis relies on the introduction of a suitable measure of 
the width of the symmetry-breaking soliton-like ground 
state characterizing the nonlinear semiclassical theory in 
the large interaction limit, which is then readily trans- 
ported to the quantum level. This quantity allows us to 
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perform a systematic comparison between the semiclas- 
sical and quantum ground state, exposing striking simi- 
larities and significantly extending the discussion in Ref. 
[lo|| . On the one hand, the comparison of the semiclassi- 
cal localized state and its quantum counterpart is made 
more quantitative. On the other hand, we extend the pa- 
rameter range to include the localization/delocalization 
transition occurring in the semiclassical nonlinear theory, 
and show that it has a clear correspondent at the quan- 
tum level. In particular, we demonstrate that the change 



in the semiclassical bifurcation pattern is maintained also 
quantum-mechanically. Our analysis is enriched by a de- 
tailed investigation of the three-site case, which makes it 
possible to visualize directly the structure of the quan- 
tum ground state. We also include a somewhat rigorous 
discussion of the relation between the quantum theory 
for attractive lattice bosons and its semiclassical version, 
related to the discrete self-trapping equations. This high- 
lights a finite-population effect. 



[1] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
[2] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[3] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and 

I. Bloch, Nature 415, 39 (2002). 
[4] L. Bernstein, Physica D 68, 174 (1993). 
[5] L. Bernstein, J. C. Eilbeck, and A. C. Scott, Nonlinearity 

3, 293 (1990). 

[6] E. Wright, J. Eilbeck, M. Hays, P. Miller, and A. Scott, 

Physica D 69, 18 (1993). 
[7] M. W. Jack and M. Yamashita, Phys. Rev. A 71, 023610 

(2005). 

[8] P. Buonsante, V. Penna, and A. Vezzani, Phys. Rev. A 

72, 043620 (2005). 
[9] N. Oelkers and J. Links, Phys. Rev. B 75, 115119 (2007). 
[10] J. Javanainen and U. Shrestha, Phys. Rev. Lett. 101, 

170405 (2008). 

[11] J. Eilbeck, P. Lomdahl, and A. Scott, Physica D 16, 318 
(1985). 

[12] A. Smerzi, A. Trombettoni, P. G. Kevrekidis, and A. R. 
Bishop, Phys. Rev. Lett. 89, 170402 (2002). 

[13] P. Buonsante, P. Kevrekidis, V. Penna, and A. Vezzani, 
J. Phys. B 39, S77 (2006). 

[14] P. Buonsante, P. G. Kevrekidis, V. Penna, and A. Vez- 
zani, Phys. Rev. E 75J316212 (2007). 

[15] The Authors of Ref. [9[] analyze the L = 3 case in de- 



tail and maintain that their results differ from those 
in Refs. [HI, [3. The differences in the observed bifur- 
cation patterns are easily explained by observing that 
what is referred to as energy in Ref. 9] is actually the 
(normalized) frequency uj appearing in the DST fixed- 
point equations J3}. The energy of the solution, as dic- 
tated by the semiclassical Hamiltonian, is given by E = 
u)N/L-\-U/2^2j |0j| 4 - Considering the frequency instead 
of the energy changes the relative ordering of some ex- 
cited states. See also Ref. [ill ]. 
[16] T.-L. Ho and C. Ciobanu, J. Low Temp. Phys. 135, 257 
(2001). 

[17] P. Ziri, J. Chwedenczuk, B. Oles, K. Sacha, and M. Trip- 

penbach, Europhysics Letters 83, 64007 (2009). 
[18] We assume unitary radius. 

[19] R. Roth and K. Burnett, PRA 68, 023604 (2003). 

[20] The static structure factor for the ring-like lattice we are 

considering would depend on a vectorial quasimomen- 

tum. The exponent in Eq. (p"2|) would be q • (r^ — r^), 

where Yj is the position of site j. 
[21] Y. Iba, Trans. Jap. Soc. for Artificial Intelligence 16, 279 

(2001). 

[22] P. Buonsante and V. Penna, J. Phys. A 41, 175301 
(2008). 

[23] K. Huang, Statistical Mechanics (John Wiley and Sons, 
1987). 



